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A single-electron tunneling (SET) device with a nanoscale central island that can move with re- 
spect to the bulk source- and drain electrodes allows for a nanoelectromechanical (NEM) coupling 
between the electrical current through the device and mechanical vibrations of the island. Although 
an electromechanical "shuttle" instability and the associated phenomenon of single-electron shuttling 
were predicted more than 15 years ago, both theoretical and experimental studies of NEM-SET struc- 
tures are still carried out. New functionalities based on quantum coherence, Coulomb correlations 
and coherent electron-spin dynamics are of particular current interest. In this article we present a 
short review of recent activities in this area. 



PACS numbers: 73.23.-b, 72.10.Fk, 73.23.Hk, 85.85.+j 
I. INTRODUCTION 

Electro-mechanical and mechano-electrical transduc- 
tions phenomena have historically contributed greatly 
to the advancement of technology in our society. To- 
day, such operations can be achieved on the single- 
molecular level with obvious advantages brought about 
by the miniaturization of the devices involved. In addi- 
tion qualitatively novel functionalities become available 
due to new physics that becomes relevant in materials 
structured on the nanometer length scale. Quantum me- 
chanics and electron-electron (Coulomb) correlations are 
defining ingredients of mesoscopic physics, which applies 
to nanoscale devices whose properties may be determined 
by a single or a few degrees of freedom. As a result there 
is a possibility and also an advantage to be gained by 
taking quantum coherence into account when designing 
electro-mechanical devices for the purpose of quantum 
manipulation and quantum communication. 

The single-electron tunneling (SET) transistor is a 
nanodevice with particularly prominent mesoscopic fea- 
tures. Here, the Coulomb blockade of single-electron tun- 
neling at low voltage bias and temperature [ ] makes 
Ohm's law for the electrical conductance invalid in the 
sense that the electrical current is not necessarily propor- 
tional to the voltage drop across the device. Instead, the 
current is due to a temporally discrete set of events where 
electrons tunnel quantum-mechanically one-by-one from 
a source to a drain electrode via a nanometer size island 
(a "quantum dot"). This is why the properties of a single 
electronic quantum state are crucial for the operation of 
the entire device. 

Since the probability for quantum mechanical tunnel- 
ing is exponentially sensitive to the tunneling distance it 



follows that the position of the quantum dot relative to 
the electrodes is crucial. On the other hand the strong 
Coulomb forces that accompany the discrete nanoscale 
charge fluctuations, which are a necessary consequence 
of a current flow through the SET device, might cause a 
significant deformation of the device and move the dot, 
hence giving rise to a strong electro-mechanical coupling. 
This unique feature makes the so-called nanoelectrome- 
chanical SET (NEM-SET) devices, where mechanical de- 
formation can be achieved along with electronic oper- 
ations, to be one of the best nanoscale realizations of 
electromechanical transduction. 

In this review we will discuss some of the latest achieve- 
ments in the nano-electromechanics of NEM-SET devices 
focusing on the new functionality that exploits quantum 
coherence in both the electronic and the mechanical sub- 
systems. The choice of materials for making a NEM-SET 
device brings an additional dimension to exploring its 
quantum performance. By choosing superconductors or 
magnets as components of the device one may, e.g., take 
advantage of a macroscopic ordering of electrons with 
respect to both their charge and spin. We will discuss 
how the electronic charge as well as the electronic spin 
contribute to electromechanical and mechano-electrical 
transduction in a NEM-SET device. New effects appear 
also due to the high mechanical deformability of molecu- 
lar NEM-SET structures (polaronic effects) which could 
be accompanied with effects of strong electron-electron 
interactions (Coulomb blockade and Luttinger liquid phe- 
nomena) as well as effects caused by strong tunneling 
coupling (Kondo- nano- mechanics). 

The phenomenon of shuttling and the sensitivity of 
electronic tunneling probabilities to mechanical deforma- 
tion of the device are in the focus of the present review. 
In this sense it is is an update of our earlier reviews of 
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shuttling [ - ]. Other aspects of nanoelectromechanics 
are only briefly discussed here. We refer readers to the 
well-known reviews of Refs. 5-9 on nanoelectromechani- 
cal systems for additional information. 

This review is organized as follows. In Section II we 
introduce the basic concepts of electron shuttling (Sub- 
section II. A) and consider the influence of polaronic ef- 
fects on the shuttle instability (Subsection II. B). In the 
end of this subsection the quantum shuttle is introduced 
and briefly discussed. Section III deals with the shuttling 
of Cooper pairs (Subsection III. A) and polaronic effects 
on the Josephson current through a vibrating quantum 
dot (Subsection III.B). A novel phenomenon — magnetic 
shuttling — is considered in Section IV, where the effects 
of spin-controlled shuttling of electric charge (Subsec- 
tion IV. A), the spintro-mechanics of the magnetic shut- 
tle (Subsection IV. B) and mechanical transportation of 
magnetization (Subsection IV. C) are discussed. In the 
end of Section III the Kondo regime of electron shuttling 
is reviewed (Subsection IV. D). Recent experiments on the 
observation of electron shuttling are briefly discussed in 
Section V. Section VI summarizes the latest theoretical 
achievements in nanoelectromechanics of shuttle devices. 
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FIG. 1: Model system consisting of a movable quantum dot 
placed between two leads. An effective elastic force acting on 
the dot due to its connections to the leads is described by a 
parabolic potential. Only one single electron state is available 
in the dot and the non-interacting electrons in the leads are 
assumed to have a constant density of states. Reprinted with 
permission from [12], D. Fedorets et al, Europhys. Lett. 58, 
99 (2002). © 2002, EDP Sciences. 



Shuttle instability in the quantum regime of 
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II. SHUTTLING OF SINGLE ELECTRONS 

A single-electron shuttle can be considered as the ul- 
timate miniaturization of a classical electric pendulum 
capable of transferring macroscopic amounts of charge 
between two metal plates. In both cases the electric 
force acting on a charged "ball" that is free to move in 
a potential well between two metal electrodes kept at 
different electrochemical potentials, eV = ijll — Hr, re- 
sults in self-oscillations of the ball. Two distinct physical 
phenomena, namely the quantum mechanical tunneling 
mechanism for charge loading (unloading) of the ball (in 
this case more properly referred to as a grain) and the 
Coulomb blockade of tunneling, distinguish the nanoelec- 
tromechanical device known as a single-electron shuttle 
[10] (see also [ ]) from its classical textbook analog. The 
regime of Coulomb blockade realized at bias voltages and 
temperatures eV, T <C Eq (where Eq — e 2 /2C is the 
charging energy, C is the grain's electrical capacitance) 
allows one to consider single electron transport through 
the grain. Electron tunneling, being extremely sensitive 
to the position of the grain relative to the bulk electrodes, 
leads to a shuttle instability — the absence of any equi- 
librium position of an initially neutral grain in the gap 
between the electrodes. 

In Subsection A below we consider the characteristic 
features of the single-electron shuttle in the case when 
the shuttle dynamics can be treated as classical motion. 
Quantum corrections to this picture are then discussed 
in Subsection B. 



Theoretically, it is convenient to study the single- 
electron shuttle in an approach [ ] where the grain 
is modeled as a single-level quantum dot (QD) that is 
weakly coupled (via a tunnel Hamiltonian) to the elec- 
trodes (see Fig. 1). The Hamiltonian corresponding to 
this model reads 

H= £ H^+H QD + £ H«\ (1) 

j=L,R j—L,R 

where the Hamiltonian 

H i j) = ^2( £ kj ~ H) a kj a kj (2) 

k 

describes noninteracting electrons in the left (j = L) and 
right (j = R) leads, which are kept at different chemical 
potential fij; a\.j(cLkj) creates (annihilates) an electron 
with momentum k in lead j. The QD Hamiltonian takes 
the form 

H QD = e cftd + SiXdU + ^(X 2 + P 2 ), (3) 

where d)(d) is the creation (destruction) operator for an 
electron on the dot, So is the energy of the resonant 
level, X is the dimensionless coordinate operator (nor- 
malized by the amplitude xq of zero-point fluctuations, 
xq = ^h/Mujo, M is the mass of QD), P is the cor- 
responding momentum operator ([X,P] = i), uq is the 
frequency of vibrons and e\ is the electromechanical in- 
teraction energy. The physical meaning of the second 
term in Eq. (3) is the interaction energy due to the cou- 
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pling of the electron charge density on the dot with the 
electric potential <p(x) = Ex crni where E = V/do is the 
electric field in the gap between electrodes (do is the dis- 
tance between the electrodes, x cm is the center-of-mass 
coordinate of the quantum dot). In this case the cou- 
pling "constant" e\ from Eq. (3) is a linear function of 
bias voltage £i(V) — eVxo/do. 

The tunneling Hamiltonian in Eq. (1) differs 

from its standard form. The explicit coordinate depen- 
dence of the tunneling matrix elements introduces ad- 
ditional electron- vibron interactions (additional to those 
described by the second term of Eq. (3). These result 

in the appearance of quantum cohesive forces = 
-dIl[ j) /dX, j = (+, -) = (L, R), where 

H t ] = Yl ^ exp(jA t X)4 + h.c. (4) 

k 

Here t$j is the bare tunneling amplitude, which corre- 
sponds to a weak dot-electrode coupling, X t = xo/h is 
a dimensionless parameter (l t is the electron tunneling 
length) that characterizes the sensitivity of the tunneling 
matrix elements to a shift of the dot center-of-mass coor- 
dinate with respect to its equilibrium (x cm = 0) position. 

In this Section we are interested in the classical shuttle 
motion (X) = x c (t) ^> 1 induced by an applied voltage 
V. Here the average (O) = Tr {p(t)0} is taken with the 
statistical operator p(t) obeying the Liouville-von Neu- 
mann equation 

iW t p{t) = [H,p{t)]. (5) 

For a weak electron- vibron interaction characterized by 
the dimensionless coupling constant A = —V2si/huo <C 
1, one can neglect the effects (~ A 2 ) of zero-point QD 
fluctuations (see below). In this case the classical dot 
coordinate x c (t) is governed by Newton's equation [12] 

x c + u%x c = F(t)/M, (6) 

where the average force F(t) = —Ti{p(t)(dH/dX)} = 
F e (t) + F c (t) consists of two terms: the electric force 
F e oc A acting on the accumulated charge on the QD, and 
the cohesive force, F c oc A £ , produced by the position- 
dependent hybridization of the electronic states of the 
grain and the leads. The evaluation of both forces can 
be done analytically [ , ] , either by solving Heisenberg 
equations of motion for the fermion operators (a^-, d) or 
by using the Keldysh Green's function approach. The 
nonlinear Eq. (6) for the classical shuttle motion can be 
analyzed in two cases: (i) near the shuttle instability 
(x c —> 0), and (ii) for the developed shuttle motion (finite 
x ci small XtX c ). 

For weak electromechanical coupling A, At <C 1 it was 
shown [ ] that the amplitude of initially small oscil- 
lations starts to grow exponentially (~ e Tst ) if eV > 
eV c = 2(^0 + ^o) 5 which means that the threshold 
voltage for the shuttle instability is V c . The rate of 
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FIG. 2: Step- like I-V curves for single-electron shuttling for 
different model parameters. Reprinted with permission from 
[12], D. Fedorets et a/., Europhys. Lett. 58, 99 (2002). © 
2002, EDP Sciences. 

growth of the instability, r s , depends on the level width 

r = r L r^/(r L + r fl ), where r, = 2^ |^| 2 K is 

the density of states, which is assumed to be an energy 
independent quantity in the wide band approximation 
[14, 15]) and on the strengths of the electromechanical 
coupling r s ~ XX t T/h. It was also shown [ ] that even 
in the absence of mechanical friction, which can be taken 
into account phenomenologically by adding to Eq. (6) 
the dissipative term jx C: the instability develops into a 
limit cycle. This is in contrast with a classical shuttle 
[10], where stability of the system can be achieved, only 
at finite mechanical dissipation 7 / 0. In the consid- 
ered model the effective dissipation is provided by the hy- 
bridization (r L (^)) of the resonant level with the metallic 
leads. 

Notice that the cohesive force F c is unimportant for 
the developed shuttle motion. However, this exchange 
interaction along with the direct electric coupling (A) de- 
termines the growth rate of the shuttle instability, which 
demonstrates the important role of electron tunneling in 
the dynamics of single electron shuttling. 

B. Strong electron-vibron interaction and 
polaronic effects in electron shuttling 

What is the role of quantum effects in the vibra- 
tional subsystem? It is known that the electron-vibron 
interaction (second term in Eq. (3)) results in vibron- 
assisted electron tunneling [ ] (the appearance of in- 
elastic channels) and for A ^> 1 in strong suppression of 
the probability of electron transfer through the elastic 
channel (Franck- Condon blockade [ ], see also the re- 
views 17, 18). Could these vibrational effects influence 
the single-electron shuttling phenomenon? 

The problem just stated was considered in Ref. 19. 
Quantum fluctuations of the dot coordinate can be taken 
into account by replacing the operator X in Eqs. (3,4) 
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FIG. 3: The increment rate of the shuttle instability as a func- 
tion of bias voltage for T/Tiujq = 0.2. The dotted line repre- 
sents the result of Ref. 3 extended to the region of strong 
electromechanical coupling. Reprinted with permission from 
[19], G.A. Skorobagat'ko et a/., Fiz. Nizk. Temp. 35, 1221 
(2009). © 2009, B. Verkin Institute for Low Temperature and 
Engineering of the NAS of Ukraine. 



by x c (t) + x, where (x) = 0. The "quantum" part of 
the electron-vibron interaction (eihx) in Eq. (3) can be 
eliminated in the Hamiltonian (4) by a standard trick 
- the Lang-Firsov unitary transformation [20, 21] U = 
exp(iApn), where [x,p] = i. After the unitary transfor- 
mation the electron-vibron coupling appears in the tun- 
neling Hamiltonian in the form of an additional opera- 
tor factor exp(zAp). Since in the transformed Hamilto- 
nian electron-vibron interactions enter only in the tun- 
neling Hamiltonian, the problem can be solved in pertur- 
bation theory with respect to the bare tunneling width 
IV In the lowest order of perturbation theory the av- 
erages of bosonic and fermionic operators are decou- 
pled and the bosonic correlation functions (exp[a*tf (i) + 
P*b(t)] exp[a&(0)+/3&t(0)])o can be evaluated analytically 
((...)o denotes an average with respect to the noninter- 
acting Hamiltonian Hq = Tiuj^b). The result [19] is an 
increment rate r s of the shuttle instability that is valid 
both for weak A <C 1 and strong A > 1 electromechanical 
coupling. For a symmetric junction Yl = Yr = Yo and 
T = it reads 



-AA,exp(-A 2 + A t 2 -AA,) ^ 
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I, 



eV 
2huj() 



1=0 



A 2 
2 



(A + A t ) 2 
11 



(7) 



(8) 



and the symbol [•] here denotes the integer part. In the 
limit of weak electromechanical couplings A <C 1, Xt <C 1 
Eq. (7) reproduces the result of Ref. , r s « YoXX t /?i, 
which predicts a linear dependence of r s on the bias volt- 
age, r s oc A oc V. Two new factors in Eq. (7) (the expo- 



nential factor and the sum over open inelastic channels) 
reflect the influence of two major vibrational effects - 
the Franck-Condon blockade and vibron-assisted electron 
tunneling - on the shuttle instability. Since the dimen- 
sionless coupling At = x$/l t is always small for molec- 
ular devices, the main quantum effect of vibrations is 
the "polaronic" narrowing of the bare tunneling width 
To — >• To exp(— X 2 (V)). It increases with the increase of 
the number of new channels (eV ^> Tiujq) and reaches its 
maximum at eVm ~ M(duo) 2 (d is the distance between 
electrodes, M is the mass of the vibrating molecule). A 
nonmonotonic dependence of the increment rate of the 
shuttle instability on V is shown in Fig. 3. One can see 
that an instability takes place in a finite interval of bias 
voltages due to the Franck-Condon blockade. Besides, at 
low temperatures, T <C hujo, the increment rate oscillates 
with a period of order Tiujq and with a relatively large am- 
plitude. Notice that at low voltages X(V) <C 1 and hence 
"polaronic" effects could not affect the "intrinsic" shuttle 
instability, which takes place at V = 2huJo/e. However, 
if the grain is pinned or if there is strong friction in the 
system, 7 ^> cjq, the shuttle instability takes place at 
much higher voltages when X(V) > 1. In this case large 
oscillations of the increment parameter could lead to un- 
usual behavior of the I-V characteristics for a shuttle- 
based single-electron transistor. A small change in bias 
voltage (smaller than hujo/e) would take the system from 
the shuttle regime of transport (with strongly enhanced 
electron tunneling probability) to the ordinary regime of 
tunnel transport (small tunneling probability) and the 
other way around. Therefore one can expect pronounced 
negative differential conductance (NDC) features (on the 
scale of hujo/e) in the current voltage characteristics (see 
also Refs. 16, 22). 

It is interesting to note that in the absence of mechan- 
ical damping, (7 = 0), the threshold voltage for single 
electron shuttling is determined by the vibron energy 
[12] eVth — 2hoJo This quantum threshold value indicates 
where the process of inelastic electron tunneling (with 
emission of a vibron) becomes energetically available. In 
the realistic case of finite friction the threshold bias volt- 
age is found by solving the equation r s (Vth) = 7, where 
r s (V) is defined in Eq. (7) (see Fig. 3). In the limit of 
weak electromechanical coupling one gets for the thresh- 
old electric field (E t h — Vth/d) the result 
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(9) 



which has a linear dependence on the rate of dissipation 
7. The exponential increase (above threshold) of the am- 
plitude A(t) of the shuttle motion means that for the 
fully developed (stationary) shuttle motion A(po) ^> x$. 
However, at the initial stage of the instability the oscilla- 
tion amplitude can be of the order of xo and the classical 
treatment of shuttle motion ceases to be valid. 

A fully quantum-mechanical approach to single elec- 
tron shuttling was developed in Refs. 23-25, where it was 
shown that the shuttle instability (exponential increase of 
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(x) and (p)) in the limit l t ^> l E = eE/Muo^ ~ #o occurs 
for a threshold electric field that coincides with Eq. (9), 
found in a quasiclassical approach. By using the Wigner 
distribution function, which allows one to visualize the 
behavior of a quantum system in phase space [ ], two 
different regimes of single electron shuttling were found. 
The classical regime (small fluctuations around the sta- 
tionary trajectory x c (t) = A c sin(co>ot)) is realized for the 
fields E > E q > E t h, where [25] 

eE q = c(^j Moj 2 l t (10) 

(C ~ 10 -2 ), i.e. in the case of weak mechanical dissi- 
pation 7 < r(xo/l t ) 4 and for large bias voltages. For 
low biases, when the electric fields, acting on the charged 
QD are in the interval E t h < E < E q , the shuttle regime 
has a specific quantum character. The Wigner function 
is strongly smeared around the classical trajectory. It is 
concentrated in a region between two circles with radii 
Rout, Rin > xo and R out - R in > x for l t > x . This 
behavior is characterized by pronounced quantum fluc- 
tuations and can be interpreted as a quantum shuttle. 
It is difficult to detect a quantum shuttle by measuring 
the average current since its qualitative behavior has no 
distinctive features in comparison with the classical shut- 
tle. The noise properties (and in general the full counting 
statistics) of NEMS will be crucial for detecting single- 
electron shuttling. 



III. MECHANICALLY MEDIATED 
SUPERCONDUCTIVITY AND POLARONIC 
EFFECTS IN THE JOSEPHSON CURRENT 

Shuttling of electric charges between nonsuperconduct- 
ing electrodes by itself does not require phase coherence. 
Even in the quantum regime of the Coulomb blockade, 
when only a single (resonant) level is involved in electron 
transport, phase coherent effects have little influence on 
electron shuttling. This is not the case for magnetic and 
superconducting leads. Magnetic exchange forces make 
the coherent electron-spin dynamics important for elec- 
tron shuttling. Superconducting transport is by defini- 
tion a phase coherent phenomenon and thus Cooper-pair 
shuttling has to be strongly different from single-electron 
shuttling. 

In this Section we consider shuttling of Cooper pairs 
between two superconducting electrodes (subsection A) 
and the influence of vibrational modes on the Josephson 
current. 



A. Shuttling of Cooper pairs 

The main requirement for the observation of a me- 
chanically mediated Josephson current [26, 27] is that 
phase coherence is preserved during the transportation 



of Cooper pairs between the two superconducting leads, 
and during the process of transferring charge between 
the bulk superconductors and a movable superconduct- 
ing grain. This requirement can be fulfilled if the super- 
conducting grain is small enough to be in the Coulomb 
blockade regime [ ] (see also [ ]) so that it can play the 
role of a single-Cooper-pair box [ ] (see also [30, 31]). 
The implication is that the characteristic energy scales 
of the small superconducting grain - the Josephson en- 
ergy Ej = (h 2 /2e)J c (where J c is the critical current) 
and the charging energy Ec — (2e) 2 /2C (where C is 
the grain capacity) - have to obey the double inequality 
Ej <C Ec <C A (where A is the superconducting gap) 
while the temperature has to be low enough to make 
T <C Ec- In this regime the single-electron states on 
the grain are energetically unfavorable (the parity effect 
[29, 32]) and the superconducting properties of the sys- 
tem can be described by a two-level model (see for exam- 
ple the review [33]). The corresponding state vector of a 
single-Cooper-pair box is a coherent superposition of the 
states with n = and n = 1 Cooper pairs ( more gen- 
erally states with the different number of Cooper pairs: 
2N and 2 (N + 1)) on the grain. For a movable Cooper 
pair box the energy scale for mechanical vibrations Huq 
has to be much smaller then all other energy scales. This 
additional requirement prevents the creation of quasipar- 
ticles and allows one to consider the mechanical motion 
of the grain as an adiabatic process. 

The Hamiltonian of the system is expressed in terms 
of the Cooper pair number operator n for the grain and 
the phases of the superconducting leads, <Pl,r : 

H = ~\ E E^x(t)}[e iVL ' R \l)(0\+h.c.}+SE c {x(t)}h. 

i=L,R 

(ii) 

The operator |1)(0| changes the number of Cooper pairs 
on the grain. An essential specific feature here is the 
dependence of the charging energy difference 5Ec = 
Ec(n = 1) — Ec(n = 0) and the coupling energies 
Ej ,R (x) = Eoexp(—SxL y R/l t ) (Sxl^r is the distance be- 
tween the grain and the respective lead) on the instanta- 
neous position x(t) of the superconducting grain. 

It is useful to separate the adiabatic motion of the 
single-Cooper-pair box between the two superconducting 
electrodes into two different parts: (i) the free motion 
( transportation region in Fig. 4), and (ii) the process 
of loading and unloading of charge near the leads ( con- 
tact region in Fig. 4). During the free motion, when the 
Josephson energy is negligibly small and the Coulomb 
term dominates, the dynamics of the qubit is reduced to 
the time evolution of the relative phase x due to the 
second term in Eq. (11), fix = 3Ec- In general the ac- 
cumulated phases are different for left-to-right and 
right-to-left motion; x± — SEct±/h. The coher- 
ent exchange of a Cooper pair between the grain and 
the lead [stage (ii)] is characterized by the dimension- 
less Josephson coupling strength 0j ~ Ejt c /Ti, where 
t c is the time spent by the grain in contact with the 
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lead (see Fig. 4). The superconducting phase difference 
V = PR — Pl, the dynamical phases x± an d the Joseph- 
son coupling strength Oj fully control the behavior of 
the Josephson current. The characteristic value of the 
mechanically assisted supercurrent (the "critical" current 
J m ) for a periodic motion of the grain (with frequency 
/ = co>o/27r) and for strong Josephson coupling Oj « 1 is 
determined by the mechanical frequency only, J m ~2ef. 
An analytical expression for the mechanically mediated 
dc Josephson current was derived in Ref. 26 and takes 
the form 



Josephson ^ 
tunneling 



Josephson 
tunneling 



J = 2ef 



sin 3 Oj cos Oj sin <£(cos <£> + cos x) 
1 — (cos 2 Oj cos x — sin 2 Oj cos <1>) 2 ' 



(12) 



where $ = cp + x+ — X- •> X = X+ + X- • The current 
Eq. (12) is an oscillating function of the superconduct- 
ing phase difference ip (see Fig. 5), which is a spectac- 
ular manifestation of a Josephson coupling between the 
remote superconductors. In the limit of weak coupling 
< 1 and vanishingly small dynamical phase (x — > 0) 
Eq. (12) is reduced to the standard Josephson formula 
J = J c sin<£, where J c ~ eEj/%. We see from Eq. (12) 
that the main qualitative effect of the dynamical phase, 
which can be controlled by the gate voltage, is a change 
of the direction of supercurrent (if cos x + cos $ < 0). For 
a given strength of the Josephson coupling the direction 
of a mechanically mediated supercurrent is determined 
by the interplay of superconducting (cp) and dynamical 
(x) phases. Notice that in Ref. it was shown that 
mechanical transportation of Cooper pairs could estab- 
lish relative phase coherence between two mesoscopic su- 
perconductors if initially they are in states with strong 
uncorrelated phase fluctuations. 



Josephson current through a vibrating quantum 
dot 



In Section 2 we considered the influence of polaronic 
effects on the electron shuttle instability. Here we an- 
alyze how vibrational degrees of freedom affect the dc 
Josephson current. When studying the transport prop- 
erties of a superconductor/quantum dot /superconductor 
(SQDS) junction we will use the same simple model for 
the QD as in the previous Section, i.e. we consider a dot 
with a single energy level So that vibrates with angular 
frequency ujq and is weakly coupled to the superconduct- 
ing leads. The latter are described by the standard BCS 
Hamiltonian with order parameter Aj = Aoe lipj (Aq is 
the superconducting gap and cpj is the phase of the or- 
der parameter for the left, j = L, and right, j = R, 
superconductor). The coupling of the dot to the leads 
is described by a tunnel Hamiltonian, which introduces 
two energy scales to the problem, viz. the partial level 
widths T^^. These are very significant for the transport 
properties. 
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FIG. 4: Illustration of the charge transport process. The cen- 
tral island moves periodically between the leads. Close to 
each turning point Cooper pair tunneling between lead and 
island is possible since the voltage on a gate electrode (not 
shown) has been set to locally remove the electrostatic en- 
ergy difference between having zero or one extra Cooper pair 
on the island (i.e., the difference in charging energy, 5Ec, is 
zero [ ]). As the island retracts from the lead, tunneling is 
exponentially suppressed (Ej = 0) while the degeneracy of 
the two charge states is lifted as the influence of the gate is 
weakened (5Ec / 0) [ ]. Reprinted with permission from 
[2], R. I. Shekhter et a/., J. Phys.: Condens. Matter 15, 
R441 (2003). © 2003, Institute of Physics and IOP Publish- 
ing Limited. 



The QD Hamiltonian reads 

H QD = €odld a + €jh(b + fr 1 ")/ V2 + Ucfi^h^ + huptfb, 

(13) 

where d cr {d\ F ) is the destruction (creation) operator for an 
electron with spin projection a =t, b n = + n^, h a = 
d\d ai b(tf) is the vibron destruction (creation) operator, 
Si is the electron- vibron interaction energy and Uc is 
the electron-electron interaction energy. In what follows 
we will assume that Ao is the largest energy scale in 
the problem. This allows one to neglect quasiparticles 
(continuum spectrum) when calculating the dc Josephson 
current. 

In general, a coupling to vibrational modes tends to 
suppress the supercurrent [35, 36]. The suppression 
mechanism is different for "hard", Tiuoq ^> T, and "soft", 
hujo <C T, vibrons. For hard vibrons and for Tiujq > Ao 
only the ground state of the vibrational subsystem is in- 
volved in Cooper pair transport through a S/QD/S junc- 
tion. Zero-point fluctuations of the QD result in strong 
(exponential) renormalization of the electron tunneling 
probability ("polaronic" narrowing of the level width). 
For strong electron- vibron coupling — i.e. for A > 1, the 
critical current is exponentially suppressed [ ], which is 
a manifestation of the Franck- Condon blockade [16, 17] 
of the supercurrent. Notice, that the Franck-Condon 
blockade will be partially removed when Tiujq <C Ao due 
to contributions of virtual side-band channels. 

The effect of the Franck-Condon blockade was 
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FIG. 5: Magnitude of the current in Eq. (12) in units of 2ef 
as a function of the phases $ and \- Regions of black corre- 
spond to no current and regions of white to J/2ef = 0.5. The 
direction of the current, is indicated in by signs (±). To best 
see the "triangular" structure of the current, the Josephson 
coupling has been chosen to be 0j = tt/3. Reprinted with 
permission from [ ], R. I. Shekhter et al., J. Phys.: Condens. 
Matter 15, R441 (2003). © 2003, Institute of Physics and 
IOP Publishing Limited. 

first predicted [35] for a nonresonant Josephson current 

J((p) = sin (/?, Jc X ^ = exp(— 2A 2 )J C (where J c is the 
critical current in the absence of electron- vibron interac- 
tion, A = 0). Later, in Ref. , it was demonstrated that 
in the case Uc = and A —> oo an analogous result, 

Jriip) = 4 X) sin(y>/2)sgn[cos(y>/2)], (14) 

4 X) = (^) exp(-A 2 ), 

holds also for the resonant (eo = 0, To = Tl = ^r) 
current. Notice the extra factor of 2 in the exponent for 
the nonresonant critical current. 

The renormalization of the bare level width is hard to 
detect in an experiment (since one would have to vary 
the strength of the electron- vibron interaction). What 
experimental manifestation of the Franck-Condon block- 
ade can one then look for? For normal transport the 
answer is that with an increase of temperature the lifting 
of the blockade is accompanied by a nonmonotonic tem- 
perature behavior of the conductance [38, 39]. An analo- 
gous behavior has been predicted [ ] for superconduct- 
ing transport, where it is the critical current that reveals 
an anomalous T-dependence. The characteristic temper- 
ature which determines the peak in a plot of J C (T) vs. T 
is determined by the polaronic energy shift E p ~ X 2 huJo- 
For T < E p the critical current increases with temper- 
ature (in the regime of temperature-enhanced Joseph- 
son coupling) while for T ^> E p the current scales as 
1/T due to a partial cancelation of Andreev levels con- 
tributions. For moderately strong electron- vibron inter- 
actions, A > 1 the crossover from low-T regime to 1/T- 



scaling looks like a "resonant" enhancement of the critical 
current at T ~ E p (see Fig. 6). 

For soft vibrons, uq — > 0, the slowly vibrating QD is 
always able to change its equilibrium position ((x) = 0) 
in order to minimize the total energy. If one neglects 
electron-electron interactions (Uc = 0) the total energy 
of the weak link can be readily evaluated in the qua- 
siclassical approximation if the dimensionless operator 
x = (b + b^)/y/2 in Eq. (13) is replaced by the classi- 
cal variable x c . Then the total energy E t consists of two 
terms: (i) the elastic energy, and (ii) the energy of the 
filled Andreev level 

Et = - ^(£0+^)2-^0082(^/2) (15) 

(for simplicity we consider here a symmetric junction). 
It is easy to see that for coupling strengths such that 
e 2 > huJoTo, a condition which is always fulfilled in the 
considered limit ujo — >• 0, the energy minimum corre- 
sponds to a shifted QD position, x c ^ 0. In an effec- 
tively asymmetric junction the resonant (eo = 0) current 
is suppressed [40] so that 

J(^) = J c sin^, Jc = ^(~^j ' ( 16 ) 

Unlike the exponential (in the electron- vibron interaction 
strength) suppression of the critical current induced by 
zero-point fluctuations of the QD coordinate (Tiuoq ^> Tq), 
soft vibrons give rise to a power-like (polaronic) suppres- 
sion of the form J oc A -2 . At finite temperatures ther- 
mally excited vibron polarons (excitations in the state 
x c 7^ 0) tend to shift the QD towards its spatially sym- 
metric position. This means that \x c (T)\ < |x c (0)|, 
which implies that at low temperatures the current grows 
with an increase of temperature. The crossover from the 
regime of a temperature enhanced supercurrent to a stan- 
dard 1/T scaling of the critical current occurs abruptly at 
T ~ T p [ ]. Therefore, both the Franck-Condon block- 
ade of the supercurrent and the polaronic effects on the 
Josephson current are manifested in an anomalous (non- 
monotonic) temperature behavior of the critical current 
(see Fig. 6). Estimations show [40] that for already exist- 
ing transport experiments on suspended single wall car- 
bon nanotubes (see e.g. [ ]) the polaronic temperature 
is in the range T p ~ (1 — 10) K, which makes the observa- 
tion of polaronic effects in carbon nanotube-based SNS 
junctions a feasible experiment. 

Finally, we discuss the influence of a charging energy 
Uc on the polaronic effects on the Josephson current. To 
this end we first note that the electron- vibron interaction 
renormalizes the electron-electron correlation energy so 
that Uc —> U^ = Uc — 2X 2 huJo and hence diminishes 
the strength of the interaction [35]. It is physically obvi- 
ous that as long as U^ < T the effects of a finite charging 
energy are negligible small. When U^ > T correlations 
split the energy level eo and for U^ ^> T the conditions 
for resonant tunneling of Cooper pairs can not be sat- 
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FIG. 6: Temperature dependence of the Josephson current 
for (p = 1, huio/si — 0.25, the value of %ujqTq/s1 — 0.2 (inset 
- huJoTo/si = 0.05). Reprinted with permission from [40], 
A. V. Parafllo et a/., (unpublished). 



isfied. If the electron-electron correlations are so strong 
that ^> Ao they additionally suppress the critical 
current by a factor T/A « 1. If E p <C Uff < A the 
charging energy does changes the value of the low-T crit- 
ical current but it can not influence the predicted anoma- 
lous temperature dependence of the current at T c± E p . 
What happenes if ~ E p is an open question, which 
needs further investigation. 



IV. ELECTRO - AND SPINTRO - MECHANICS 
OF MAGNETIC SHUTTLE DEVICES 

In this Section we will explore new functionalities that 
emerge when nanomechanical devices are partly or com- 
pletely made of magnetic materials. The possibility of 
magnetic ordering brings new degrees of freedom into 
play in addition to the electronic and mechanical ones 
considered so far, opening up an exciting perspective to- 
wards utilising magneto-electro-mechanical transduction 
for a large variety of applications. Device dimensions in 
the nanometer range mean that a number of mesoscopic 
phenomena in the electronic, magnetic and mechanical 
subsystems can be used for quantum coherent manipu- 
lations. In comparison with the electromechanics of the 
nanodevices considered above the prominent role of the 
electronic spin in addition to the electric charge should 
be taken into account. 

The ability to manipulate and control spins via electri- 
cal [42-44] magnetic [ ] and optical [46] means has gen- 
erated numerous applications in metrology [47] in recent 
years. A promising alternative method for spin manip- 
ulation employs a mechanical resonator coupled to the 
magnetic dipole moment of the spin(s), a method which 
could enable scalable quantum information architectures 



[48] and sensitive nanoscale magnetometry [49-51]. Mag- 
netic resonance force microscopy (MRFM) was suggested 
as a means to improve spin detection to the level of 
a single spin and thus enable three dimensional imag- 
ing of macromolecules with atomic resolution. In this 
technique a single spin, driven by a resonant microwave 
magnetic field interacts with a ferromagnetic particle. 
If the ferromagnetic particle is attached to a cantilever 
tip, the spin changes the cantilever vibration parameters 
[52]. The possibility to detect [52] and monitor the co- 
herent dynamics of a single spin mechanically [ ] has 
been demonstrated experimentally. Several theoretical 
suggestions concerning the possibility to test single-spin 
dynamics through an electronic transport measurement 
were made recently [54-57]. Complementary studies of 
the mechanics of a resonator coupled to spin degrees of 
freedom by detecting the spin dynamics and relaxation 
were suggested in [54-61] and carried out in [ ]. Elec- 
tronic spin-orbit interaction in suspended nanowires was 
shown to be an efficient tool for detection and cooling 
of bending-mode nanovibrations as well as for manipu- 
lation of spin qubit and mechanical quantum vibrations 
[63-65]. 

An obvious modification of the nano-electro-mechanics 
of magnetic shuttle devices originates from the spin- 
splitting of electronic energy levels, which results in the 
known phenomenon of spin-dependent tunneling. Spin- 
controlled nano-electro-mechanics which originates from 
spin-controlled transport of electric charge in magnetic 
NEM systems is represented by number of new magneto- 
electro-mechanical phenomena. 

Qualitatively new opportunities appear when magnetic 
nanomechanical devices are used. They have to do with 
the effect of the short-ranged magnetic exchange inter- 
action between the spin of electrons and magnetic parts 
of the device. In this case the spin of the electron rather 
than its electrical charge can be the main source of the 
mechanical force acting on movable parts of the device. 
This leads to new physics compared with the usual elec- 
tromechanics of non-magnetic devices, for which we use 
the term spintro- mechanics. In particular it becomes pos- 
sible for a movable central island to shuttle magnetiza- 
tion between two magnetic leads even without any charge 
transport between the leads. The result of such a me- 
chanical transportation of magnetization is a magnetic 
coupling between nanomagnets with a strength and sign 
that are mechanically tunable. 

In this Section we will review some early results that 
involve the phenomena mentioned above. These only 
amount to a first step in the exploration of new oppor- 
tunities caused by the interrelation between charge, spin 
and mechanics on a nanometer length scale. 



A. Spin-controlled shuttling of electric charge 

By manipulating the interaction between the spin of 
electrons and external magnetic fields and/or the inter- 
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nal interaction in magnetic materials, spin-controlled na- 
noelectromechanics may be achieved. 

A new functional principle — spin-dependent shuttling 
of electrons — for low magnetic field sensing purposes 
was proposed by Gorelik et al. in Ref. . This principle 
may lead to a giant magnetoresistance effect in external 
magnetic fields as low as 1-10 Oe in a magnetic shuttle de- 
vice if magnets with highly spin-polarized electrons (half 
metals [67-71]) are used as leads in a magnetic shuttle de- 
vice. The key idea is to use the external magnetic field to 
manipulate the spin of shuttled electrons rather than the 
magnetization of the leads. Since the electron spends a 
relatively long time on the shuttle, where it is decoupled 
from the magnetic environment, even a weak magnetic 
can rotate its spin by a significant angle. Such a rota- 
tion allows the spin of an electron that has been loaded 
onto the shuttle from a spin-polarized source electrode 
to be reoriented in order to allow the electron finally to 
tunnel from the shuttle to the (differently) spin-polarized 
drain lead. In this way the shuttle serves as a very sen- 
sitive "magnetoresistor" device. The model employed in 
Ref. assumes that the source and drain are fully po- 
larized in opposite directions. A mechanically movable 
quantum dot (described by a time-dependent displace- 
ment x(t)), where a single energy level is available for 
electrons, performs driven harmonic oscillations between 
the leads. The external magnetic field, i7, is perpendicu- 
lar to the orientations of the magnetization in both leads 
and to the direction of the mechanical motion. 

The spin-dependent part of the Hamiltonian is speci- 
fied as 

H magn (t) = J{t){a\a^-a\a^)-^^{a\a^a\a^)^ (17) 

where J(t) = Jr(£) — JlW? Jl{r){^) are the ex- 
change interactions between the on-grain electron and 
the left (right) lead, g is the gyromagnetic ratio and /i is 
the Bohr magneton. The proper Liouville-von Neumann 
equation for the density matrix is analyzed and an av- 
erage electrical current is calculated for the case of large 
bias voltage. 

In the limit of weak exchange interaction, J max <C \±H 
one may neglect the influence of the magnetic leads on 
the on-dot electron spin dynamics. The resulting current 
is 

eujp sin 2 (^/2) tanh(w/4) 
7T sin 2 (tf/2) +tanh 2 (w/4) 

where w is the total tunneling probability during the con- 
tact time to, while $ ~ irgfiH/hujQ is the rotation angle 
of the spin during the "free-motion" time. 

The theory [66] predicts oscillations in the magnetore- 
sistance of the magnetic shuttle device with a period 
AH p , which is determined from the equation Tuujq = 
<7/x(l + w)AH p . The physical meaning of this relation 
is simple: every time when ujq/Q = n + 1/2 (Q = gjiH/Ti 
is the spin precession frequency in a magnetic field) the 



shuttled electron is able to flip fully its spin to remove 
the "spin-blockade" of tunneling between spin polarized 
leads having their magnetization in opposite directions. 
This effect can be used for measuring the mechanical fre- 
quency thus providing dc spectroscopy of nanomechani- 
cal vibrations. 

Spin-dependent shuttling of electrons as discussed 
above is a property of non-interacting electrons, in the 
sense that tunneling of different electrons into (and out 
of) the dot are independent events. The Coulomb block- 
ade phenomenon adds a strong correlation of tunneling 
events, preventing fluctuations in the occupation of elec- 
tronic states on the dot. This effect crucially changes 
the physics of spin-dependent tunneling in a magnetic 
NEM device. One of the remarkable consequences is 
the Coulomb promotion of spin-dependent tunneling pre- 
dicted in Ref. . In this work a strong voltage depen- 
dence of the spin-flip relaxation rate on a quantum dot 
was demonstrated. Such relaxation, being very sensi- 
tive to the occupation of spin-up and spin-down states 
on the dot, can be controlled by the Coulomb block- 
ade phenomenon. It was shown in Ref. 72 that by lift- 
ing the Coulomb blockade one stimulates occupation of 
both spin-up and spin-down states thus suppressing spin- 
flip relaxation on the dot. In magnetic devices with 
highly spin-polarized electrons electronic spin-flip can be 
the only mechanism providing charge transport between 
oppositely magnetized leads. In this case the onset of 
Coulomb blockade, by increasing the spin-flip relaxation 
rate, stimulates charge transport through a magnetic 
SET device (Coulomb promotion of spin-dependent tun- 
neling). Spin-flip relaxation qualitatively also modifies 
the noise characteristics of spin-dependent single-electron 
transport. In Refs. r 3, 74 it was shown that the low- 
frequency shot noise in such structures diverges as the 
spin relaxation rate goes to zero. This effect provides an 
efficient tool for spectroscopy of extremely slow spin-flip 
relaxation in quantum dots. Mechanical transportation 
of a spin-polarized dot in a magnetic shuttle device pro- 
vides new opportunities for studying spin-flip relaxation 
in quantum dots. The reason can be traced to a spin- 
blockade of the mechanically aided shuttle current that 
occurs in devices with highly polarized and colinearly 
magnetized leads. As was shown in Ref. the above 
effect results in giant peaks in the shot-noise spectral 
function, wherein the peak heights are only limited by 
the rates of electronic spin flips. This enables a nanome- 
chanical spectroscopy of rare spin-flip events, allowing 
spin-flip relaxation times as long as 10 /is to be detected. 

The spin-dependence of electronic tunneling in mag- 
netic NEM devices permits an external magnetic field 
to be used for manipulating not only electric transport 
but also the mechanical performance of the device. This 
was demonstrated in Refs. 76, 77. A theory of the quan- 
tum coherent dynamics of mechanical vibrations, electron 
charge and spin was formulated and the possibility to 
trigger a shuttle instability by a relatively weak magnetic 
field was demonstrated. It was shown that the strength 



of the magnetic field required to control nanomechanical 
vibrations decreases with an increasing tunnel resistance 
of the device and can be as low as 10 Oe for giga-ohm 
tunnel structures. 

A new type of nanoelectromechanical self excitation 
caused entirely by the spin splitting of electronic en- 
ergy levels in an external magnetic field was predicted 
in Ref. 78 for a suspended nanowire, where mechanical 
motion in a magnetic field induces an electromotive cou- 
pling between electronic and vibrational degrees of free- 
dom. It was shown that a strong correlation between 
the occupancy of the spin-split electronic energy levels in 
the nanowire and the velocity of flexural nanowire vibra- 
tions provides energy supply from the source of DC cur- 
rent, flowing through the wire, to the mechanical vibra- 
tions thus making possible stable, self-supporting bend- 
ing vibrations. Estimations made in Ref. show that 
in a realistic case the vibration amplitude of a suspended 
carbon nanotube (CNT) of the order of 10 nm can be 
achieved if magnetic field of 10 T is applied. 



B. Spintro- mechanics of magnetic shuttle device 

New phenomena, qualitatively different from the elec- 
tromechanics of nonmagnetic shuttle systems, may ap- 
pear in magnetic shuttle devices in a situation when 
short-range magnetic exchange forces become compara- 
ble in strength to the long-range electrostatic forces be- 
tween the charged elements of the device [ ]. There is 
convincing evidence that the exchange field can be several 
tesla at a distance of a few nanometers from the surface 
of a ferromagnet [1 -82] . Because of the exponential de- 
cay of the field this means that the force experienced by a 
single-electron spin in the vicinity of magnetic electrodes 
can be very large. These spin-dependent exchange forces 
can lead to various "spintro-mechanical" phenomena. 

Mechanical effects produced by a long-range electro- 
static force and short-ranged exchange forces on a mov- 
able quantum dot are illustrated in Fig. 7. The elec- 
trostatic force acting on the dot, placed in the vicinity 
of a charged electrode (Fig. 7(a)), is determined by the 
electric charge accumulated on the dot. In contrast, the 
exchange force induced by a neighboring magnet depends 
on the net spin accumulated on the dot. While the elec- 
trostatic force changes its direction if the electric charge 
on the dot changes its sign, the spin-dependent exchange 
force is insensitive to the electric charge but it changes 
direction if the electronic spin projection changes its sign. 
A very important difference between the two forces is that 
the electrostatic force changes only as a result of injec- 
tion of additional electrons into (out of) the dot while the 
spintronic force can be changed due to the electron spin 
dynamics even for a fixed number of electrons on the dot 
(as is the case if the dot and the leads are insulators). In 
this case interesting opportunities arise from the possibil- 
ity of transducing the dynamical variations of electronic 
spin (induced, e.g., by magnetic or microwave fields) to 
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FIG. 7: A movable quantum dot in a magnetic shuttle de- 
vice can be displaced in response to two types of force: (a) 
a long-range electrostatic force causing an electromechanical 
response if the dot has a net charge, and (b) a short-rang mag- 
netic exchange force leading to "spintromechanical" response 
if the dot has a net magnetization (spin). The direction of 
the force and displacements depends on the relative signs of 
the charge and magnetization, respectively. Reprinted with 
permission from [83], R. I. Shekhter et aL, Phys. Rev. B 86, 
100404 (2012). © 2012, American Physical Society. 



mechanical displacements in the NEM device. In Ref. 83 
a particular spintromechanical effect was discussed - a gi- 
ant spin-filtering of the electron current (flowing through 
the device) induced by the formation of what we shall 
call a "spin-polaronic state". 

The Hamiltonian that describes the magnetic nanome- 
chanical SET device in Ref. 83 has the standard form (its 
spin-dependent part depends now on the mechanical dis- 
placement of the dot). Hence H = Hi ead +H tunne i+H dot , 
where Hi eada = J2k 7 *,s e ksaa{ S(T aksa describes electrons 
(labeled by wave vector k and spin a =t,l) in the two 
leads (s = L,R). Electron tunneling between the leads 
and the dot is modeled as 

Htunnel = ^ T ^ X ) a ksa C ^ + H ' C ' ( 19 ) 
k,a,s 

where the matrix elements T s (x) = exp(=fx/l t ) (k 
is the characteristic tunneling length) depend on the dot 
position x. The Hamiltonian of the movable single-level 
dot is 

H dot = hu b t b^^2n (T [eo-sgn(a)J(x)}^Ucri t n i , (20) 

cr 

where sgn(t, X) = ±1, Uc is the Coulomb energy associ- 
ated with double occupancy of the dot and the eigenval- 
ues of the electron number operators n a is or 1. The 
position dependent magnitude J(x) of the spin dependent 
shift of the electronic energy level on the dot is due to the 
exchange interaction with the magnetic leads. Here we 
expand J(x) to linear order in x so that J(x) = +jx 
and without loss of generality assume that = 0. 

The modification of the exchange force, caused by 
changing the spin accumulated on the dot, shifts the equi- 
librium position of the dot with respect to the magnetic 
leads of the device. Since the electron tunneling ma- 
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Spatial position of the dot 

FIG. 8: Diagram showing how the equilibrium position of 
the movable dot depends on its net charge and spin. The 
difference in spatial displacements discriminates transport 
through a singly occupied dot with respect to the electron 
spin. Reprinted with permission from [ ], R. I. Shekhter et 
al, Phys. Rev. B 86, 100404 (2012). © 2012, American 
Physical Society. 



trix element is exponentially sensitive to the position of 
the dot with respect to the source and drain electrodes 
one expects a strong spin-dependent renormalization of 
the tunneling probability, which exponentially discrim- 
inates between the contributions to the total electrical 
current from electrons with different spins. This spatial 
separation of dots with opposite spins is illustrated in 
Fig. 8. While changing the population of spin-up and 
spin-down levels on the dot (by changing e.g. the bias 
voltage applied to the device) one shifts the spatial posi- 
tion x of the dot with respect to the source/drain leads. 
It is important that the Coulomb blockade phenomenon 
prevents simultaneous population of both spin states. If 
the Coulomb blockade is lifted the two spin states be- 
come equally populated with a zero net spin on the dot, 
S = 0. This removes the spin-polaronic deformation and 
the dot is situated at the same place as a non-populated 
one. In calculations a strong modification of the vibra- 
tional states of the dot, which has to do with a shift of its 
equilibrium position, should be taken into account. This 
results in a so-called Frank- Condon blockade of electronic 
tunneling [16, 17]. The spintro-mechanical stimulation of 
a spin-polarized current and the spin-polaronic Franck- 
Condon blockade of electronic tunneling are in compe- 
tition and their interplay determines a non-monotonic 
voltage dependence of the giant spin-filtering effect. 

To understand the above effects in more detail con- 
sider the analytical results of Ref. . A solution of the 
problem can be obtained by the standard sequential tun- 
neling approximation and by solving a Liouville equation 
for the density matrix for both the electronic and vibronic 
subsystems. The spin-up and spin-down currents can be 
expressed in terms of transition rates (energy broaden- 
ing of the level) and the occupation probabilities for the 
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FIG. 9: Spin polarization of the current through the model 
NEM-SET device under discussion. Reprinted with permis- 
sion from [ ], R. I. Shekhter et a/., Phys. Rev. B 86, 100404 
(2012). © 2012, American Physical Society. 



dot electronic states. For simplicity we consider the case 
of a strongly asymmetric tunneling device. At low bias 
voltage and low temperature the partial spin current is 



exp 



_ ( x o y 



sgn(a)/3 , (21) 



where j3 = x^/Tiuoolt. In the high bias voltage (or tem- 
perature) regime, max{eV 1 T ^> E p }, where the pola- 
ronic blockade is lifted (but double occupancy of the dot 
is still prevented by the Coulomb blockade), the current 
expression takes the form 



h 



exp [2n B + 1] 



■2sgn(a)0 



(22) 



where ub is Bose-Einstein distribution function. The 
scale of the polaronic spin-filtering of the device is de- 
termined by the ratio /3 of the polaronic shift of the 
equilibrium spatial position of a spin-polarized dot and 
the electronic tunneling length. For typical values of the 
exchange interaction and mechanical properties of sus- 
pended carbon nanotubes this parameter is about 1-10. 
As was shown this is enough for the spin filtering of the 
electrical current through the device to be nearly 100 % 
efficient. The temperature and voltage dependence of the 
spin- filtering effect is presented in Fig. 9. The spin filter- 
ing effect and the Franck-Condon blockade both occur at 
low voltages and temperatures (on the scale of the pola- 
ronic energy; see Fig. 9 (a)). An increase of the voltage 
applied to the device lifts the Franck-Condon blockade, 
which results in an exponential increase of both the cur- 
rent and the spin- filtering efficiency of the device. This 
increase is blocked abruptly at voltages for which the 
Coulomb blockade is lifted. At this point a double occu- 
pation of the dot results in spin cancellation and removal 
of the spin-polaronic segregation. This leads to an ex- 
ponential drop of both the total current and the spin 
polarization of the tunnel current (Fig. 9 (b)). As one 
can see in Fig. 9 prominent spin filtering can be achieved 
for realistic device parameters. The temperature of oper- 
ation of the spin-filtering device is restricted from above 
by the Coulomb blockade energy. One may, however, 
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FIG. 10: Single- domain magnetic grains with magnetic mo- 
ments Ml and M# are coupled via a magnetic cluster with 
magnetic moment m, the latter being separated from the 
grains by insulating layers. The gate electrodes induce an 
ac electric field, concentrated in the insulating regions. This 
field, by controlling the heights of the tunnel barriers, affects 
the exchange magnetic coupling between different components 
of the system. Reprinted with permission from [ ], L. Y. 
Gorelik et al., Phys. Rev. Lett. 91, 088301 (2003). © 2003, 
American Physical Society. 



consider using functionalized nanotubes [ ] or graphene 
ribbons [ ] with one or more nanometer-sized metal or 
semiconductor nanocrystal attached. This may provide 
a Coulomb blockade energy up to a few hundred kelvin, 
making spin filtering a high temperature effect [83]. 



vectors in such a way that the vector sum is conserved. 
This is why the result of this rotation can be viewed as a 
transfer of some magnetization Am from one ferromag- 
net to the other. As a result the magnetization of the dot 
experiences some rotation around a certain axis. The to- 
tal angle <fi of the rotation accumulated during the time 
when the dot is magnetically coupled to the lead is an es- 
sential parameter which depends on the mechanical and 
magnetic characteristics of the device. The continuation 
of the mechanical motion breaks the magnetic coupling 
of the dot with the first lead but later, as the dot ap- 
proaches the other magnetic lead an exchange coupling 
is established with this second lead with the result that 
magnetization which is "loaded" on the dot from the first 
lead is "transferred" to the this second lead. This is how 
the transfer of magnetization from one magnetic lead to 
another is induced mechanically. The transfer creates 
an effective coupling between the magnetizations of the 
two leads. Such a non-equilibrium coupling can be effi- 
ciently tuned by controlling the mechanics of the shuttle 
device. It is particularly interesting that the sign of the 
resulting magnetic interaction is determined by the sign 
of cos(0/2). Therefore, the mechanically mediated mag- 
netic interaction can be changed from ferromagnetic to 
anti-ferromagnetic by changing the amplitude and the 
frequency of mechanical vibrations [86]. 



C. Mechanically assisted magnetic coupling 
between nanomagnets 

The mechanical force caused by the exchange interac- 
tion represents only one effect of the coupling of magnetic 
and mechanical degrees of freedom in magnetic nano- 
electromechanical device. A complementary effect is the 
of mechanical transportation of magnetization, which we 
are going to discuss in this subsection. 

In the magnetic shuttle device presented in Fig. 10, 
a ferromagnetic dot with total magnetic moment m is 
able to move between two magnetic leads, which have 
total magnetization M^^. Such a device was suggested 
in Ref. 86 in order to consider the magnetic coupling 
between the leads (which in their turn can be small mag- 
nets or nanomagnets) produced by a ferromagnetic shut- 
tle. It is worth to point out that the phenomenon we 
are going to discuss here has nothing to do with trans- 
ferring electric charge in the device and it is valid also 
for a device made of nonconducting material. The main 
effect, which will be in the focus of our attention, is the 
exchange interaction between the ferromagnetic shuttle 
(dot) and the magnetic leads. This interaction decays 
exponentially when the dot moves away from a lead and 
hence it is only important when the dot is close to one of 
the leads. During the periodic back-and-forth motion of 
the dot this happens during short time intervals near the 
turning points of the mechanical motion. An exchange 
interaction between the magnetizations of the dot and 
a lead results in a rotation of these two magnetization 



D. Resonance spin-scattering effects. Spin shuttle 
as a "mobile quantum impurity". 

The Kondo effect in electron tunneling results from the 
spin exchange between electrons in the leads and the is- 
land (quantum dot) that couples the leads and manifests 
itself as a sharp zero bias anomaly in the low-temperature 
tunneling conductance. Many-particle interactions and 
the tunneling renormalize the electron spectrum enabling 
Kondo resonances both for odd [87] and even [88, 89] elec- 
tron occupations. In the latter case the Kondo resonance 
is caused by the singlet-triplet crossover in the ground 
state (see [{ ] for a review). In the simplest case of odd 
occupancy a cartoon of a quantum well and a schematic 
Density of States (DoS) is shown in Fig. 11. For simplic- 
ity we consider a case when the dot is occupied by one 
electron (as in a SET transistor). The dot level is not 
in resonance with the Fermi level of the leads (ep), but 
located at an energy —Ed, below it. The dot is in the 
Coulomb blockade regime and the corresponding charg- 
ing energy is denoted as Ec- The resonance spin scat- 
tering results in the formation of a narrow peak in the 
DoS known as the Abrikosov-Suhl resonance [91-93] (see 
Fig. 11, right panel). The width of this resonance defines 
a unique energy scale, the Kondo temperature TV, which 
determines all thermodynamic and transport properties 
of the SET device through a one-parametric scaling [93] . 
The width Y of the dot level, associated with the tun- 
neling of dot electrons to the continuum of levels in the 
leads, is assumed to be smaller than the charging energy 
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FIG. 11: Nanomechanical resonator with spin as a "mobile 
quantum impurity". 



Ec, providing a condition for an integer valency regime. 
When the shuttle moves between source (S) and drain 
(D) (see the lower panel of Fig. 11), both the energy Ed 
and the width T acquire a time dependence. This time 
dependence results in a coupling between mechanical and 
electronic degrees of freedom. If a source-drain voltage 
V s d is small enough (eV s d <C Tk) the charge degree of 
freedom of the shuttle is frozen out while spin plays a 
very important role in co-tunneling processes. Namely, 
the dot electron's spin can be flipped while the electron 
tunnels from the left to the right lead. Thus, the initial 
and final states of the quantum impurity can have differ- 
ent spins. This process is accompanied by simultaneous 
creation of spin excitations in the Fermi sea. The many- 
electron scattering processes then lead to the formation 
of an Abrikosov-Suhl resonance. This resonance can be 
viewed as a Kondo cloud built up from both conduction 
electrons in the leads and a localized electron in the dot. 
Since all electrons in the cloud contain information about 
the same impurity, they are mutually correlated. Thus, 
NEM providing a coupling between mechanical and elec- 
tronic degrees of freedom introduces a powerful tool for 
manipulation and control of the Kondo cloud and gives 
a very promising and efficient mechanism for electrome- 
chanical transduction on the nanometer length scale. 

Building on an analogy with the shuttling experiments 
of Refs. and , let us consider a device where an 
isolated nanomachined island oscillates between two elec- 
trodes (Fig. 11, lower panel). The applied voltage is as- 
sumed low enough so that the field emission of many 
electrons, which was the main mechanism of tunneling 
in those experiments, can be neglected. We emphasize 
that the characteristic de Broglie wave length associated 
with the dot should be much shorter than typical dis- 
placements allowing thus for a classical treatment of the 
mechanical motion of the nano-particle. The condition 
hujQ <C Tk, necessary to eliminate decoherence effects, 
requires for e.g. planar quantum dots with the Kondo 



temperature Tk ^ 100 mK, the condition uoq < 1 GHz 
for oscillation frequencies to hold; this frequency range is 
experimentally feasible [94, 95]. The shuttling island is 
then to be considered as a "mobile quantum impurity", 
and transport experiments will detect the influence of 
mechanical motion on the differential conductance. If the 
dot is small enough, then the Coulomb blockade guaran- 
tees the single electron tunneling or cotunneling regime, 
which is necessary for the realization of the Kondo effect 
[30, 96]. Cotunneling is accompanied by a change of spin 
projection in the process of charging/discharging of the 
shuttle and therefore is closely related to the spin/charge 
pumping problem [97]. 

A generic Hamiltonian for describing the resonance 
spin-scattering effects is given by the Anderson model, 

#o = Xl^ aC L,a c ^a + y^i E d ~ eEx]dl a d i(T + E c n 2 

k,a i<J 

H tun = J2 T i i} ( x ) \<L, Jia + H.c] , (23) 

ika,a 

where c\ a , d\ a create an electron in the lead a=L,R, or 
the dot level £^=1,2, respectively, n = d} i(T d i(T , E is the 
electric field between the leads. The tunnelling matrix el- 
ement Tj^ R {x) = T^'ft exp[=Fx(t)/Zt], depends exponen- 
tially on the ratio of the time-dependent displacement 
x(t) and the electronic tunnelling length l t . The time- 
dependent Kondo Hamiltonian can be obtained from it 
by applying a time-dependant Schrieffer- Wolff transfor- 
mation [98]: 



H 



K 



koto ,k' a,' <j' 



Jaa>(t)[a aa 'S + -C'jc^a^'' 



',a' (24) 



where J a ^(t) = y/T a (t)T af (t)/(7rp E d (t)) and S = 
\d\G a(7 id! a ^ T a (t) = 27rpo\T a (x(t))\ 2 are level widths due 
to tunneling to the left and right leads. 

As long as the nano-particle is not subject to an exter- 
nal time-dependent electric field, the Kondo temperature 
is given by T K = Dq exp [— (ttEc)/(STo)] (for simplicity 
we assumed that Tl(0) = T^(0) = To; Dq plays the role 
of effective bandwidth). As the nano-particle moves adi- 
abatically, Tiujq <C To, the decoherence effects are small 
provided Tiujq <C T k . 

Let us first assume a temperature regime T ^> Tk 
(weak coupling). In this case we can build a per- 
turbation theory controlled by the small parameter 
p$J(t) Ih[Dq/T] < 1 assuming time as an external pa- 
rameter. The series of perturbation theory can be 
summed up by means of a renormalization group proce- 
dure [ ]. As a result, the Kondo temperature becomes 
oscillating in time: 



T K (t)=D(t)exp 



7tE C 



(25) 



8r cosh(2x(t)// t )_ 
Neglecting the weak time-dependence of the effective 
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FIG. 12: Differential conductance G of a Kondo shuttle 
for which Yo/U=0A. The solid line denotes G for a shut- 
tle with Tl=Fr, A=lt, the dashed line shows G for a static 
nano- island with Yl = I\r, A=0, the dotted line gives G for 
rL/rn=0.5, A=0. The inset shows the temporal oscillations 
(here Q = cjo) of Tk for small A=0. 05 Zt (dotted line) and 
large A=2.bl t (solid line) shuttling amplitudes. Reprinted 
with permission from [99], M. N. Kiselev et al, Phys. Rev. B 
74, 233403 (2006). © 2006, American Physical Society. 
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FIG. 13: Time dependence of the current Io for different 
values of asymmetry parameter u = xo/lt. Here red, blue 
and black curves correspond to u — 0.5;1.0;1.5;. For all 
three curves shuttle oscillates with amplitude x m ax — h, 
tuo /(T K ) min = 10" 3 , |eVb ias |/TF n = 9VBB/TZ in = 0.1 
with T { ° } = 2K, h/L = 10" 4 . Reprinted with permission 
from [ ], M. N. Kiselev et al Phys. Rev. Lett. 110, 
066804 (2013). © 2013, American Physical Society. 



bandwidth D(t) ~ Do, we arrive at the following ex- 
pression for the time- averaged Kondo temperature: 



(Tk) 



rriO 
1 K 



exp 



ttE c smh 2 (x(t)/l t ) 



_4T l + 2smh 2 (x(t)/l t )_ 



. (26) 



Here (...) denotes averaging over the period of the 
mechanical oscillation. The expression (26) acquires 
an especially transparent form when the amplitude of 
the mechanical vibrations A is small: A < l t . In 
this case the Kondo temperature can be written as 
(Tk) = T£exp(-2jy), with the Debye- Waller-like ex- 
ponent W = —7rEc(x 2 (t)))/(8r lt), giving rise to the 
enhancement of the static Kondo temperature. 

The zero bias anomaly (ZBA) in the tunneling conduc- 
tance is given by 



G(T) = ^-Go( 



4r L (t)r fl (t) 



1 



(r L (t) + r*(t))2 [HT/T K (t)W 



where Go = e 2 /h is a unitary conductance. Although 
the central position of the island is most favorable for the 
Breit-Wigner (BW) resonance (Tl = T^), it corresponds 
to the minimal width of the Abrikosov-Suhl resonance. 
The turning points correspond to the maximum of the 
Kondo temperature given by the equation (25) while the 
system is away from the BW resonance. These two com- 
peting effects lead to the effective enhancement of G at 
high temperatures (see Fig. 12). 

Summarizing, it was shown in [ ] that Kondo shut- 
tling in a NEM-SET device increases the Kondo temper- 
ature due to the asymmetry of coupling at the turning 
points compared to at the central position of the island. 
As a result, the enhancement of the differential conduc- 
tance in the weak coupling regime can be interpreted as 
a pre-cursor of strong electron-electron correlations ap- 



pearing due to formation of the Kondo cloud. 

Let us consider now the strong coupling regime, T <C 
Tk- The current through the system subject to a con- 
stant source-drain bias V s d can be separated in two parts: 
a dc current associated with a time-dependent dc conduc- 
tance and an ac current related to the periodic motion 
of the shuttle. For both currents Kondo physics plays an 
important role. While the dc current is mostly responsi- 
ble for the frequency shift, the ac current gives an access 
to the dynamics of the Kondo cloud and provides infor- 
mation about the kinetics of its formation. In order to 
evaluate both contributions to the total current we ro- 
tate the electronic states in the leads in such a way that 
only one combination of the wave functions is coupled 
to the quantum impurity. The cotunneling Hamiltonian 
may be rationalized by means of the Glazman-Raikh ro- 
tation, paramet rized by the a ngle d t defined by the rela- 
tion tantf, = ^\Y R (t)/Y L (t)\. 

Both the ac and dc contributions to the current can 
be calculated by using Noziere's Fermi-liquid theory (see 
[ ] for details). The ac contribution, associated with 
the time dependence of the Friedel phase 5 a [ ] , is given 
by 



Lc(t) = 



x(t) eE c eV sd 



tanh 



"77 8fr7 ' T K (t) ' cosh 2 ^2[^-xol^ 



(28) 



(exp(4x //t) = r i? (0)/r L (0)) and the "ohmic" dc contri- 
bution is fully defined by the adiabatic time-dependence 
of the Glazman-Raikh angle 

I D c(t) = G V sd sin 2 2$ t sin 2 S a (29) 

a 

As a result, the ac contribution to the total current can 
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be considered as a first non-adiabatic correction: 



hot = hd(x(t)) - x 



. dl ad HttEc 



16r T^ 0) 



(30) 



where I ad = 2 • G • ^ cosh- 2 (2[x(t) - ar ]/Z t ) and T^ 0) 
is the Kondo temperature at the equilibrium position. 
The small correction to the adiabatic current in (30) 
may be considered as a first term in the expansion over 
the small non adiabatic parameter uot <C 1, where r 
is the retardation time associated with the inertia of 
the Kondo cloud. Using such an interpretation one gets 
r = h7rE c /(WT TP). 

Equation (30) allows one to obtain information about 
the dynamics of the Kondo clouds from an analysis of an 
experimental investigation of the mechanical vibrations. 
The retardation time associated with the dynamics of the 
Kondo cloud is parametrically large compared with the 
time of formation of the Kondo cloud tk = Ti/Tk and 
can be measured owing to a small deviation from adia- 
baticity. Also we would like to emphasize a supersensi- 
tivity of the quality factor to a change of the equilibrium 
position of the shuttle characterized by the parameter u 
(see Fig. 13). The influence of strong coupling between 
mechanical and electronic degrees of freedom on the me- 
chanical quality factor has been considered in [101]. It 
has been shown that both suppression Q > Qo and en- 
hancement Q < Qo of the dissipation of nanomechanical 
vibrations (depending on external parameters and the 
equilibrium position of the shuttle) can be stimulated by 
Kondo tunneling. The latter case demonstrates the po- 
tential for a Kondo induced electromechanical instability. 

Summarizing, we emphasize that the Kondo phe- 
nomenon in single electron tunneling gives a very 
promising and efficient mechanism for electromechani- 
cal transduction on a nanometer length scale. Measur- 
ing the nanomechanical response on Kondo-transport in 
a nanomechanical single-electron device enables one to 
study the kinetics of the formation of Kondo-screening 
and offers a new approach for studying nonequilibrium 
Kondo phenomena. The Kondo effect provides a possi- 
bility for super high tunability of the mechanical dissi- 
pation as well as super sensitive detection of mechanical 
displacement. 



V. EXPERIMENTAL OBSERVATION OF 
ELECTRON SHUTTLING 

Electron shuttling was theoretically predicted to occur 
in mechanically soft mesoscopic systems about 15 year 
ago [ ]. Since then there has been a steadily increasing 
interest in studying this nonequilibrium electromechani- 
cal phenomenon from both theoretical and experimental 
points of view. Each year the technical capability to fab- 
ricate shuttle-like devices improves. On the experimental 
side there are two main directions in the study of mechan- 
ically mediated electron transport: (i) electron shuttling 




FIG. 14: Electron micrograph of a "quantum bell": The pen- 
dulum is clamped on the upper side of the structure. It can be 
set into motion by an ac-power, which is applied to the gates 
on the left- and right-hand side (Gl and G2) of the clapper 
(C). Electron transport is then observed from source (S) to 
drain (D) through the island on top of the clapper. The island 
is electrically isolated from the rest of the clapper, which is 
grounded. Reprinted with permission from [102], A. Erbe et 
a/., Appl. Phys. Lett. 73, 3751 (1998). © 1998, American 
Institute of Physics. 



in NEMS with intrinsic electromechanical coupling, and 
(ii) electron transfer caused by an external excitation 
of mechanical motion. Here we briefly review several re- 
cent publications, which have claimed to observe electron 
shuttling. We start with the nanoelectromechanical de- 
vices based on vibrating cantilevers. 

Externally driven nanomechanical shuttles have been 
designed in Refs. 95, 102. In these experiments a nanome- 
chanical pendulum was fabricated on a Si-on-insulator 
substrate using electron and optical lithography. A metal 
island was placed on a clapper, which could vibrate be- 
tween source and drain electrodes (see Fig. 14). The 
pendulum was excited by applying an ac voltage between 
two gates on the left- and right-hand sides of the clap- 
per. The observed tunneling source-drain current was 
strongly dependent on the frequency of the exciting sig- 
nal having pronounced maxima at the eigenfrequencies 
of the mechanical modes. This fact signalizes a shuttling 
mechanism of electron transfer at typical shuttle frequen- 
cies of about 100 MHz. The measured average dc current 
at 4.2 K corresponds to 0.11 ± 0.001 electrons per cycle 
of mechanical motion. Both a theoretical analysis and 
numerical simulations showed that a large portion of the 
voltage also acts on the island. 

A very important modification of the setup in Fig. 14 
was presented in Ref. 103. There a silicon cantilever 
is part of a mechanical system of coupled resonators, 
which is a construction that makes it possible to drive 
the shuttle mechanically with a minimal destructive in- 
fluence from the actuation dynamics on the shuttle itself. 
This is achieved by a clever design that minimizes the 
electrical coupling between the driving part of the device 
(either a magnetomotively driven, doubly clamped beam 
resonator, or a capacitively coupled remote cantilever) 
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FIG. 15: (a) SEM micrograph and experimental circuitry of 
the silicon nanopillar studied in Ref. : At the source (S), 
an ac signal, V ac , is applied with a superimposed dc bias 
V. The net current Id is detected at the drain (D) with 
a current amplifier. The third electrode (G) is floating, (b) 
Finite-element simulation of the base oscillation mode that 
compiles for the nanopillar to fo = 5367 MHz. (c) When 
the island is deflected toward one electrode, the instantaneous 
voltage bias determines the preferred tunneling direction. Co- 
tunneling is absent in this case due to an increased distance to 
the opposite electrode. Reprinted with permission from [105], 
D. V. Scheibe and R. H. Blick, Appl. Phys. Lett 84, 4632 
(2004). © 2004, American Institute of Physics. 



and the driven part (the cantilever that carries the shut- 
tle on its tip). In principle, systems of this type can be 
used for studies of shuttle transport through supercon- 
ducting and magnetic systems. 

In Ref. 104 the role of nanocantilever was played by a 
semiconductor nanowire. In was shown that under cer- 
tain conditions the constant electron beam produced by 
a scanning electron microscope can excite self-sustained 
mechanical oscillations of semiconducting SiC nanowire. 
The nanowire plays the role of mechanical resonator and 
may be represented by an RC circuit element (where R is 
the nanowire resistance and C is the capacitance between 
the nanowire end and its environment or some electrode 
placed near the resonator). The periodic electrostatic 
force, which depends on the charge of the wire, acts on 
the wire due to variations in capacitance. The charge 
on the nanowire is also a time-dependent function. Dis- 
charging occurs due to RC relaxation accompanied by 
the drift of electrons to the tungsten tip, on which the 
wire is attached. The charging is provided by the electron 
beam (for details see [ ]). The semiconductor nanowire 
starts to oscillate and goes to the stationary cycle. Thus 
in [104] a new type of electromechanical coupling was 
studied. 

Interesting results on mechanically assisted charge 
transfer were obtained in Ref. 105 for a device fabricated 
as a silicon nanopillar located between source and drain 
contacts (see Fig. 15). The device was manufactured in a 




FIG. 16: Two coupled electron shuttles realized as nanopil- 
lars. The metallic top layer allows electron exchange with the 
source (S) and drain (D) contacts. The scale bar corresponds 
to a length of 200 nm. The inset shows a broader view of the 
coplanar- waveguide into which the nanopillars are embedded. 
The scale bar in the inset is 10 /mi. Reprinted with permis- 
sion from [ ], C. Kim et a/., ACS Nano 6, 651 (2012). © 
2012, American Chemical Society. 



two step process: nanoscale lithography using a scanning 
electron microscope (SEM) and, second, dry etching in 
a fluorine reactive ion etcher (RIE). The lithographically 
defined gold structure acts both as electrical current leads 
and etch mask for the RIE. A simple geometry defined 
by SEM consequently results in the freestanding isolat- 
ing nanopillar of intrinsic silicon with a conducting metal 
(Au) island at its top (see Fig. 15). This island serves as 
the charge shuttle. The metal island and the nanopillar 
are placed in the center of two facing electrodes. The 
system is biased by an ac voltage at source, rather than 
a sole dc bias, to avoid the dc-self excitation. Application 
of an ac-signal excites one of the nanopillar eigenmodes 
resonantly. The device was operated at room tempera- 
ture and the capacitance was not sufficiently small to re- 
alize the Coulomb blockade regime. The dependencies of 
the current on bias frequency, as well as on an additional 
dc bias, allowing to tune resonances, were measured. 

In Refs. 106 and 107 electron transport at room tem- 
perature through two nanopillars was considered. The 
nanopillars act as shuttles placed in series between source 
and drain electrodes under ac/dc excitation (see Fig. 16). 
The linear size of the island on top of each pillar was 
65 nm with pillar heights of 250 nm and an inter-pillar 
distance of 17 nm. At first the I-V characteristics at room 
temperature were measured [ ] without any ac signal. 
The bias voltage dependence of the current was shown to 
be almost linear (except for a small deviation ex V 2 above 
1 V). At low voltages (< 100 mV) Coulomb blockade fea- 
tures (a CB staircase) are apparent. From experimental 
data one finds the charging energy of the two coupled 
nanopillars to be Eq = 41 meV, which is larger than 
the (room) temperature equivalent of the experimental 
setup. 

The dependence of current on gate and bias voltages 
is manifested in the Coulomb diamonds (see Fig. 17) in 
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FIG. 17: Experimental Coulomb blockade (CB) diamonds 
traced in the normalized current in color representation. The 
lower borders of the CB regions are represented in red with 
the Coulomb plateaus depicted in green. The borders of the 
CB determined from the theory plots are marked by dashed 
black lines as a guide to the eye. Reprinted with permission 
from [ ], C. Kim et a/., ACS Nano 6, 651 (2012). © 2012, 
American Chemical Society. 
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FIG. 18: Full frequency sweep of the direct current through 
the coupled shuttle revealing the mechanical mode structure 
for Vdc = 0. Reprinted with permission from [106], C. Kim et 
al, Phys. Rev. Lett. 105, 067204 (2010). © 2010, American 
Physical Society. 



agreement with theoretical calculations. Here the typical 
Coulomb blockade staircase (superimposed on an ohmic 
response) is smeared due to thermal broadening and shut- 
tling effects. 

Then, a radio frequency signal (1 MHz - 1 GHz) was 
fed without any dc bias added. Due to the alternating- 
voltage induced ac current one would expect a zero- 
average (over an ac period) dc current without the me- 
chanical subsystem However, in the considered experi- 
ment a nonzero net current was observed. The authors 
explained this fact by the excitation of a mechanical mo- 
tion of the nanopillars. The effect of current rectification 
indicates a dynamical violation of P-symmetry (symme- 
try with respect to coordinates reflection r — »• — r) in 
the system. The direction and the amplitude of the 
dc current depends on the ac frequency (see Fig. 18). 
The presence of a broad set of resonances indicates the 




FIG. 19: I-V traces for two mechanical modes of the shuttle 
device shown in Fig. 16. As can be seen the slopes at the 
plateaus (indicated by dashed lines) increase from a shuttling 
frequency of 285 MHz (red squares) to 500 MHz (blue circles). 
Reprinted with permission from [ ], C. Kim et a/., ACS 
Nano 6, 651 (2012). © 2012, American Chemical Society. 



existence of different mechanical modes of the coupled 
nanopillars (Fig. 18). Different mechanical modes were 
investigated by applying a low dc bias voltage. The ob- 
served I-V characteristics (with the zero-bias current sub- 
tracted) is plotted in Fig. 20. The step-like features in the 
current- voltage dependence can be interpreted [ ] as a 
signature of electron shuttling in the Coulomb blockade 
regime. 

Another trend in the study of electron shuttles is to 
mimic shuttling effects by time dependent tunneling bar- 
riers. The authors of Ref. 108 investigated single-electron 
"shuttling" with a silicon metal-oxide-semiconductor 
quantum dot at low temperature (300 mK). The ana- 
lyzed system represents an electron layer at a Si/Si02 
interface below an aluminium top gate. This layer forms 
a conduction channel between source and drain electrodes 
(see Fig. 19). The quantum dot is formed by two addi- 
tional gates, which produce tunnel barriers. The effects 
of mechanical degrees of freedom of the quantum dots are 
mimicked here by an ac voltage of frequency f p applied 
to both gates. Periodic variations of the barrier heights 
induce periodic effective "displacements" of the quantum 
dot on the nanometer scale. Fig. 21 demonstrates the de- 
pendence of the current on the bias voltage for different 
frequencies, Id c = ^nef p . The observed results can be 
explained by a sequential tunneling model with the elec- 
tron temperature is a fitting parameter. Electrons are 
"hotter" than the environment due to ac voltage heating. 

An elegant experiment was performed by Park et al. 
[109], who studied electron transport through a fullerene 
molecule placed in the gap between two gold electrodes. 
The experimental results (a stair-case-like dependence of 
current on bias voltage) are well explained by the process 
of vibron-assisted tunneling [ )] . The shuttle model was 
also used for an explanation of the step-like features on 
I-V characteristics (see the discussion in reviews [2, 3]). 
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FIG. 20: (a) Schematic cross section of a fabricated silicon 
metal-oxide-semiconductor (MOS) quantum dot [ ]. Two 
aluminum barrier gates (BL and BR) are below a top gate 
(TOP) isolated with A\ x O y . The source and drain are ther- 
mally diffused with phosphorus, (b) Scanning electron mi- 
croscope image of the device with a simplified measurement 
setup, (c) Energy landscape through the dot and lead reser- 
voirs with an illustration of the electron shuttling. When 
the sinusoidal ac voltage on BR, Vbr achieves its maximum 
(dashed green line), an electron tunnels into the dot from the 
right electron reservoir. After one half of an operation pe- 
riod Vbl is at its maximum value (dotted purple line) and 
the electron tunnels away. Blue regions denotes the states in 
the leads occupied by electrons. Reprinted with permission 
from [ ], K. W. Chan et a/., Appl. Phys. Lett. 98, 212103 
(2011). © 2009, American Institute of Physics. 
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FIG. 22: Experimental (symbols) and simulated (dotted and 
dashed lines) current- volt age characteristics for shuttle junc- 
tions fabricated A. V. Moskalenko et al. [ ]. The dotted 
line corresponds to oscillations in the case of zero pinning and 
the dashed lines to the case of finite pinning in the system. 
Insets show the shuttle displacement as a function of time 
for two points, one of which is below and the other is above 
the transition into the shuttling regime. The leakage current 
through a monolayer of octanedithiol molecules is shown by 
the dashed-dotted line. Reprinted with permission from [111], 
A. V. Moskalenko et a/., Phys. Rev. B 79, 241403 (2009). © 
2009, American Physical Society. 
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FIG. 23: Confirmation that current flows through the 
nanoparticle in the shuttling regime (curve 1). If the nanopar- 
ticle is removed from the gap using an AFM tip the result is 
a drop in the current through the device of several orders of 
magnitude (curve 2). The inset shows hysteretic behavior of 
I-V curves obtained for a working shuttle junction in regimes 
of increasing and decreasing applied voltage. Reprinted with 
permission from [111], A. V. Moskalenko et a/., Phys. Rev. B 
79, 241403 (2009). © 2009, American Physical Society. 



FIG. 21: (a) Charge shuttling current measured for different 
frequencies, (b) Measured charge shuttling current (solid blue 
line) for f p = 60 MHz. The measured current is compared 
with a simulation (dotted red line) based on a sequential tun- 
neling model with variable tunneling resistances in the bar- 
riers. Reprinted with permission from [108], K. W. Chan et 
a/., Appl. Phys. Lett. 98, 212103 (2011). © 2009, American 
Institute of Physics. 



In Ref. Ill the authors reported shuttling by a 20 nm 
gold particle pasted into the gap (10-20 nm) between two 
electrodes and attached to them through a monolayer of 
organic molecules (1.8 octanedithiol). Experimental I- 
V characteristics are presented in Fig. 22. One can see 
that the theoretical fit using shuttle model is in a good 
agreement with experiment. At low voltages the cur- 
rent through the device is absent (the nanoparticle is in 



19 




FIG. 24: A nanomechanical electron shuttle, (a) SEM im- 
age of a shuttle indicating the dimensions of the gold island 
(yellow) and the high stress silicon nitride resonator (red) sus- 
pended above the silicon substrate (grey) . Tilted view of gold 
island between source and drain. Reprinted with permission 
from [ ],D. R. Koenig and E. M. Weig, Appl. Phys. Lett 
101, 213111 (2012).© 2009, American Institute of Physics. 



a locked state). At high voltages the nanoparticle be- 
gins to vibrate (the quantum dot escapes from the locked 
state) , because more electrons are transferred to the gran- 
ule and the electrostatic force acting on it becomes strong 
enough to de-pin the shuttle. The onset of shuttle insta- 
bility occurs at higher bias voltages than predicted by 
a theory assuming frictionless shuttle motion. This dis- 
crepancy can be explained by the binding of the particle 
to the electrodes (respectively the value of the threshold 
voltage is increased). Experimental data in favor of this 
assumption is the observation of an I-V hysteresis loop 
(see Fig. 23). Without a gold particle the current through 
the device is determined by the sequential tunneling of 
electrons. Comparison of the experimental data and the- 
oretical calculations suggests that the average number of 
electrons that are involved in the shuttling at voltages 
of order 3 V is about 20. This experiment can be inter- 
preted as electron shuttling but the device is operated in 
a regime very far from single electron tunneling (Coulomb 
blockade regime). Is it possible to fabricate a shuttle-like 
device operated in the Coulomb blockade regime by dc 
bias voltage? Although this shuttle was not observed yet, 
the "road" to its discovery now is sufficiently clear. 

In Ref. 112 the authors claimed to have found a 
charge shuttle (by the observation of sustained self- 
oscillation) operated solely by an applied dc voltage with- 
out external actuation. The island of the shuttle device 
in Ref. is a gold particle with typical dimensions 
35 nmx270 nmx40 nm, placed at the center of a sus- 
pended silicon nitride string (see Fig. 24). Experiments 
were performed in helium exchange gas with a pressure of 
0.5 mbar in a helium dewar at T = 4.2 K. Forty- four sam- 
ples were studied for statistical processing. To observe 
voltage-sustained self-oscillations the resonant acoustic 
drive is turned off when the island is charged. After that 
the stable charge transport is observed for almost 2000 s 
(about 10 10 shuttle cycles). The reasons for the collapse 
of self-sustained oscillations may be different: due to 
impact-induced coupling to out-of-plane or torsional mo- 
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FIG. 25: DC voltage-sustained electron shuttling and back- 
ground current. (a) Current- volt age curves of voltage- 
sustained self-sustained oscillation. Both the blue and the 
red trace, corresponding to downward and upward voltage 
sweeps, respectively, feature a sharp dissipation threshold, 
(b) Background current determined by measuring Isd in the 
absence of mechanical shuttling as a function of bias volt- 
age. The dashed box indicates the voltage range depicted in 
(a), where the background current is also shown as a black 
line. Reprinted with permission from [112], D. R. Koenig and 
E. M. Weig, Appl. Phys. Lett 101, 213111 (2012).© 2009, 
American Institute of Physics. 



tion or wear-induced alternation of island and electrodes, 
etc. In the absence of mechanical shuttling the current 
is a linear function of bias voltage (see Fig. 25). The 
presence of electron shuttling is indicated by a step-like 
current dependence. The average number of transferred 
electrons (N c± 200) is easily found from the I-V char- 
acteristics (see Fig. 25). The advantage of the dc-biased 
self-sustained shuttle current is the significant drop of 
external heat load on the system. The described experi- 
ment opens the pathway to observe a Coulomb blockade 
shuttle operated solely by a dc bias. 

One of the most promising objects to be used for ob- 
serving electron shuttling is a suspended carbon nan- 
otube. To observe an electron shuttle it is important 
to fabricate an electromechanical system with two prop- 
erties: (i) a soft mode of mechanical vibrations that is 
susceptible to becoming unstable when a bias voltage is 
applied, and (ii) a sufficiently strong electromechanical 
coupling that can overcome mechanical damping. In de- 
vices based on suspended carbon nanotubes both these 
conditions are realized [41, 113]. Three different vibra- 
tional modes (the radial mode uo = const, the stretching 
mode ujq oc 1/L and the bending mode ujq oc 1/L 2 , L 
is the nanotube length) can be excited in transport ex- 
periments with suspended carbon nanotubes. Recently it 
was shown that in a suspended carbon nanotube actuated 
by rf-radiation, the mechanical subsystem is sensitive to 
single-electron transport [ ] (see Fig. 26). In Ref. 113 
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Single electron tuning. 


(A) 


Nanotube 



versus gate voltage showing single-electron tunneling at the 
peaks and Coulomb blockade in the valleys. (B) Normal- 
ized resonance signal versus rf frequency and gate voltage 
(V s d = 1.5 mV). The tuned mechanical resonance shows up as 
the darker curve with dips at the Coulomb peaks. Reprinted 
with permission from [113], G. A. Steele et a/., Science 325, 
1103 (2009). © 2009, American Association for the Advance- 
ment of Science. 



the influence of a gate voltage on the frequency of vibra- 
tional modes in the Coulomb blockade regime of electron 
transport was studied. The bending mode was excited 
by radiation of a nearby rf antenna at the resonance fre- 
quency. This frequency was shown to drop at bias volt- 
ages such that the Coulomb blockade is lifted and an 
additional electron is transferred through the nanotube. 
The experiment demonstrates strong electromechanical 
coupling in suspended carbon nanotubes. 



VI. CONCLUSIONS 

This review demonstrates that during the last several 
years there has been significant activity in the study of 
nanoelectromechanical (NEM) shuttle structures. New 
physics was harvested for theoretical suggestions of how 
to achieve new functionality of shuttle devices. A short 
list of these suggestions would include the following 
items: 

1) The polaronic approach for studying nanoelectrome- 
chanics beyond the limit of weak NEM coupling was de- 
veloped. This approach, resulting in a qualitative modi- 
fication of the voltage and temperature control on NEM 
performance, could be applied to molecular-based shut- 
tle devices known for their high mechanical deforma- 
bility. 2) New NEM phenomena are possible in super- 
conducting shuttle devices. The superconducting cur- 
rent through Josephson weak links supported by Andreev 
bound states [ ] (see also [115], [116]) is qualitatively 
modified by coupling the latter to nanomechanical vibra- 
tions. Both mechanically assisted Cooper-pair transfer 
and electronically assisted cooling of a nanomechanical 
subsystem become possible. Strong electron-vibron in- 
teractions are manifested in the Franck- Condon blockade 
of the critical current (T = 0) and its nonmonotonous 
(anomalous) temperature behavior. 3) The possibility 



to achieve an external control on the direction of energy 
transfer between electronic and mechanical subsystems 
has been demonstrated. Efficient ground state cooling 
of a nanomechanical resonator was predicted in both su- 
perconducting and non-superconducting NEM-SET sys- 
tems [117, 118]. 4) The role of the electronic spin is 
important in shuttle devices made of magnetic materials. 
Spin-dependent exchange forces can be responsible for a 
qualitatively new nanomechanical performance opening a 
new field of study that can be called spintro- mechanics. 
5) Electronic many-body effects, appearing beyond the 
weak tunneling approach, result in single electron shut- 
tling assisted by Kondo-resonance electronic states. The 
possibility to achieve a high sensitivity to coordinate dis- 
placement in electromechanical transduction along with 
the possibility to study the kinetics of the formation of 
many-body Kondo states has also been demonstrated re- 
cently. 

There are still a number of unexplored shuttling 
regimes and systems, which one could focus on in the 
nearest future. In addition to magnetic and supercon- 
ducting shuttle devices one could explore hybrid struc- 
tures where the source/drain and gate electrodes are hy- 
brids of magnetic and superconducting materials. Then 
one could expect spintromechanical actions of a super- 
current flow as well as superconducting proximity effects 
in the spin dynamics in magnetic NEM devices. An ad- 
ditional direction is the study of shuttle operation un- 
der microwave radiation. In this respect microwave as- 
sisted spintromechanics is of special interest due to the 
possibility of microwave radiation to resonantly flip elec- 
tronic spins. As in ballistic point contacts [116] such 
flips can be confined to particular locations by the choice 
of microwave frequency, allowing for external tuning of 
the spintromechanical dynamics of the shuttle. Polaronic 
effects in superconducting and magnetic shuttle devices 
represents another interesting direction for the future re- 
search. Promising nano-objects where charge shuttling 
could play a significant role include bio-systems. Elec- 
tromechanical phenomena involving both electrons and 
protons could be very important for "transport" effects 
in cells and bacteria (see [119-121]). 

Experimental studies of shuttle devices have been de- 
veloped significantly during the last couple of years. A 
number of new implementations of the idea of nanome- 
chanical shuttling has been suggested (see chapter V). 
The aim of having fully controllable mechanics of the 
shuttle motion in combination with the lowest possible 
level of dissipation motivates the modifications which 
have been made. The possibility to achieve mechan- 
ically controllable shuttling of a single electron was a 
challenge for a number of years. Evidence indicating 
that a Coulomb blockade shuttle device has been made 
was reported recently [ ]. Although the regime of self 
supported shuttle vibrations has been reached in a num- 
ber of experiments (see Section 5), the observation of a 
shuttle instability still remains a challenge for the future. 
The achievements made in the technology of producing 
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shuttle devices [112] give promises for the construction of 
magnetic and superconducting shuttle structures, which 
would enable one to explore a number of tempting func- 
tionalities theoretically suggested during recent years. 
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